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Abstract. We use evolutionary calculations presented in 
a recent paper (Cassisi et al. 1997a: hereinafter Paper I) 
to predict B,V,I magnitudes for stars in galactic globulars. 
\ The effect of the adopted mixing length on stellar mag- 
. nitudes and colors is discussed, showing that the uncer- 
I tainty on such a theoretical parameter prevents the use of 
i MS stars as bona fide theoretical standard candles. How- 
I ever, comparison with Hipparcos data for field subdwarfs 
■ discloses a substantial agreement between theory and ob- 
I servation. Present predictions concerning the magnitude 
. of TO and of HB stars are compared with similar results 
] appeared in the recent literature. We find that our pre- 
r dictions about the dependence on metallicity of ZAHB 
magnitudes appear in good agreement with observational 
. constraints as recently discussed by Gratton et al. (1997c). 
' We present and discuss a theoretical calibration of the dif- 
. ference in magnitude between HB and TO as evaluated 
' with or without element sedimentation. The effect of a 
. variation of the original helium content on the magnitude 
I of MS, TO and HB stars is explored and discussed. Finally 
1 we use theoretical HB magnitudes to best fit the CM di- 
agram of M68 and M5, taken as representative of metal 
poor and intermediate metallicity galactic globulars, de- 
riving an age of llibl.O Gyr and 10±1.0 Gyr, respectively, 
for the adopted chemical compositions, plus an additional 
uncertainty of ±1.4 Gyr if the uncertainty on the chemical 
composition is taken into account. This result is discussed 
on the basis of current evaluations concerning cluster ages 
and distance moduli. 
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1. Introduction 

Following the effort undertaken by Hubble Space Tele- 
scope for improving the determination of the Hubble con- 
stant and by taking advantage of several improvements 
of the input physics a number of authors have recently 
revisited the prediction of stellar evolutionary theory to 
re-evaluate the age of galactic globular clusters (see, e.g., 
VandenBerg et al. 1996; D Antona et al. 1997; Salaris et al. 
1997; Brocato et al. 1997; Chaboyer et al. 1998; Salaris & 
Weiss 1998a, b). In a previous paper (Cassisi et al. 1998: 
hereinafter Paper I) we discussed on purely theoretical 
grounds the effect of the improved physical inputs on stel- 
lar evolutionary models, evaluating the additional influ- 
ence of element diffusion on the evolutionary history of 
such stars. In this paper we take advantage of such an 
homogeneous set of stellar models to present theoretical 
predictions concerning observational magnitudes and col- 
ors, and to discuss the theoretical scenario leading to the 
evaluation of cluster ages. 

Such age evaluations, as primarily based on the cali- 
bration of the Turn-Off luminosity in terms of the cluster 
age, requires the adoption of suitable standard candles 
belonging to the cluster, as given either by MS or by HB 
stars. In the last case one is dealing with the so called AV 
method, which is only formally independent of the clus- 
ter distance modulus. In all cases, connecting theory to 
ages requires several steps worth being discussed in some 
details. Thus, presenting the evolutionary results, we will 
also investigate the degree of reliability concerning theo- 
retical predictions about the magnitude of MS, TO and 
HB stars. The main aim of this investigation is to point 
out theoretical uncertainties but also to show how reli- 
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able the new models are, testing the results with available 
observational constraints and, in particular, with the ab- 
solute magnitudes of metal poor stars recently provided 
by Hipparcos. 

2. Prom theory to observation 

In order to cover with sufficient details the range in metal- 
licity Z=2-I0~'' - 6-I0^'^ which appears adequate for galac- 
tic globulars, computations presented in Paper I have been 
supplemented with an additional set of stellar models with 
Z=0.0006 . As in Paper I, evolutionary models have been 
computed by adopting an amount of original He Y=0.23 
and with a mixing length parameter 1 = 1.6 Hp. Data for 
these new models are reported in the Appendix of this 
paper. 

On the basis of these new results, in this paper we will 
rely on a grid of evolutionary models covering, at least, 
the range of cluster ages 8 -18 Gyr for the assumed metal- 
licities Z= 0.0002, 0.0006, 0.001, 0.006. The computations 
have been performed both by including or neglecting el- 
ement sedimentation. In addition, for Z=0.001 we com- 
puted two set of evolutionary tracks with element sedi- 
mentation, under the same assumption about the mixing 
length but with Y= 0.21 or 0.25, respectively. These mod- 
els will allow to discuss the influence of original He on 
models where element diffusion is taken into account. Ac- 
cording to the usual procedure, the evolutionary models 
have been used to derive cluster isochrones in the logL, 
logTe plane. 

To allow a comparison with observation, this theoret- 
ical scenario has to be implemented with suitable pro- 
cedures transforming the theoretical quantities L (lumi- 
nosity) and Te (effective temperature) into reliable pre- 
dictions about magnitudes in selected photometric bands. 
To look into this problem Fig. |l| shows our Z=0.0002 12 
Gyr isochrone (with diffusion) as transformed into the CM 
diagram according to selected available choices about the 
color-temperature and magnitude luminosity transforma- 
tions. One finds that the most recent evaluations given by 
Castelli et al. (1997 a,b) are very similar to previous eval- 
uations given by Buser & Kurucz (1978, 1992) all along 
the TO region as well as in the redder portion of the HB, 
whereas Kurucz (1992) gives slightly redder TO colors. In 
all cases one finds that the difference affects mainly the 
colors, whereas bolometric corrections appear remarkably 
similar. Thus V magnitudes appears independent of the 
choice of the models. 



Fig. 1. The run of the theoretical 12 Gyr isochrone in the 
V,B-V diagram for the labeled choices about the adopted 
model atmospheres. 



If not otherwise advised, in this paper we will rely on 
the model atmospheres by Castelli et al. (1997 a,b). Here 
we note that preliminary computations (Castelli, private 
communication) for the case Z=0.0002 have already shown 
that the model atmospheres can be further improved. As 
a matter of fact, one finds that when both the new solar 
abundances and the enhancement of a-elements is taken 
into account, model atmospheres for metal poor stars 
([Fe/H]= -2.0) reach a better agreement with available 
empirical estimate by Alonso et al. (1996). Fortunately, 
the same Fig. |l| shows that such improvement affects only 
the slope of the MS at the lower luminosity, leaving in 
particular unaffected the predicted magnitude of candi- 
date candle stars. 

As is well known, MS models are also affected by an 
intrinsic uncertainty due to the assumptions about the 
value of the mixing length parameter. The effect of varying 
the adopted mixing length within the (reasonable) interval 
1= 1.3-2.0 Hp is shown in Fig. ^, where we report selected 
isochrones for the labeled values of ages and mixing length 
parameters. One finds that for ages ranging from 9 to 15 
Gyr the MS color at Mv=6.0 decreases by about A(B- 
V)~0.02 mag, whereas decreasing the mixing length from 
2.0 down to 1.3 Hp (the local pressure scale height) gives 
again a decrease of this color of about A(B-V)~0.03 mag, 
depending on the cluster age. Thus the color of the MS at 
Mv=6.0 is far from being a firm prediction of theory. The 
same figure shows that the magnitude of the MS at B- 
V=0.55 (i.e. around theoretical predictions for Mv=6) for 
a given cluster metallicity can move by about 0.28 mag.(~ 
0.18 mag. for a variation from 1.3 to 2.0 Hp and ~ 0.1 mag. 
for a difference in age from 9 to 15 Gyr), preventing the use 
of theoretical MS models as precise distance calibrators. 

When moving toward more advanced evolutionary 
phases, as already noticed by Chaboyer (1995), the cal- 
ibration of TO magnitudes in terms of cluster ages does 
depend on the assumptions about the mixing length pa- 
rameter. One finds that such a dependence increases when 
the cluster ages and/or the metallicity is increased, for the 
simple reason that in both cases TO stars become cooler 
and, thus, more affected by external convection. As a re- 
sult, the most metal poor and younger TO (i.e., the hottest 
TO) are barely affected by the treatment of convection all 
over the explored range of mixing length parameters (1.3 
2.0 Hp). As an example, for Z=0.0002 and for the age 
t= 9 Gyr, a change from 1= 1.3 to 2.0 Hp moves the TO 
magnitude by only 0.02 mag. However, for the same as- 
sumed metallicity, at t=14 Gyr, the same variation of the 
mixing length parameter gives a difference in the TO mag- 
nitude of the order of 0.1 mag. 

According to the quoted evidences, we will proceed in 
the discussion of theoretical predictions with the warn- 
ing that discussed magnitudes should be reliable, whereas 
colors are affected by the discussed uncertainties in both 
stellar and atmospheric models. Detailed data for the 
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Fig. 2. The 9, 12 and 15 Gyr isochrones as computed under 
the labeled assumptions about the mixing length parameter 
and translated in the CM diagram according to Castelli et al. 
(1997 a,b). 

computed isochrones are available by anonymous ftp at 
astr 1 8pi . difi. unipi . it (/pub/globular ) . 

In the next section we discuss low mass stars during the 
major phase of H burning, namely from the initial main se- 
quence to the subgiant phase, just after the exhaustion of 
central H. We selected this evolutionary phase according 
to the evidence that the corresponding stellar models rely 
on the "minimum" input physics as given by an equation 
of state (EOS) not too far from a perfect gas, by radiative 
opacities and, finally, by sufficiently well known H burn- 
ing rates (see, e.g., Brocato et al. 1998) Correspondingly, 
the quoted evolutionary phases have to be regarded as the 
most solid predictions of theory. Section 4 will deal with 
a discussion of more advanced evolutionary phases, from 
the red giant branch to the AGB through the phase of cen- 
tral He burning, where more physics has to be added, like 
neutrino production by weak interactions or the physical 
behavior of electron-degenerate matter. In both cases we 
find that theoretical predictions reach a good agreement 
with observational constraints. Following such an agree- 
ment, we discuss relevant evolutionary features and, in 
particular, the estimate of cluster ages based on the HE 
luminosity level. A brief final discussion closes the paper. 

3. H burning phases. 

As a first test of theoretical results. Fig. || (a to c) shows a 
comparison between isochrones (with sedimentation) and 
Hipparcos data (Gratton et al. 1997, Chaboyer et al. 1998) 
for field subdwarfs in selected intervals of [Fe/H] values. 
Each panel compares isochrones at t=10 and 12 Gyr and 
for a given [Fe/H] value with subdwarfs with [Fe/H] esti- 
mates within ±0.1 dex. Triangles in each panel show the 
estimated shift in the isochrone if the suggested enhance- 
ment in a-elements [Q;/Fe]~0.3 (Carretta & Gratton 1997, 
Gratton et al. 1998a) is accounted as an increase in the 
total value of Z. For the sake of comparison, each panel 
shows also the location of the MS loci but for different 
metallicities. 

Inspection of the three panels shows that a-elements 
play a minor role, since the corresponding shift of the 
isochrone lies within the observational errors. Bearing in 
mind that sedimentation does not affect the location of 
MS stars, one finds that the agreement between theory 
and observation appears as fine as it can be expected ac- 
cording to the standard errors in absolute magnitudes and 
the spread in metallicity. One finds that the popular Reid's 
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statement " current models overestimate the change in lu- 
minosity with decreasing abundance for extreme metal 
poor subdwarf is hardly supported. In fact, the three 
panels of Fig.^ suggest that no contradiction appears be- 
tween extant theory and Hipparcos observations. 



Fig. 3. Comparison of present isochrones with Hipparcos sub- 
dwarf absolute magnitudes (Gratton et al. 1997). Each panel 
shows the run of the 10 and 12 Gyr isochrones (full lines) for 
the given value of [Fe/H]. Triangles (and the dashed line) show 
the shift of the 10 Gyr isochrones caused by an enhancement 
[a/Fe] ~ 0.3. As a comparison, dotted lines show the run of 10 
Gyr isochrones with different [Fe/H], as labeled. Filled squares 
and open squares indicate single stars and detected or sus- 
pected binary stars, respectively. 

Figure ^ shows the comparison between present pre- 
dictions of B-V colors at Mv= 6.0 and Hipparcos obser- 
vational data as presented by Gratton et al. (1997). One 
finds that theoretical data appear in satisfactory agree- 
ment with observations, whose error appears of the same 
order of magnitude of the theoretical uncertainty dis- 
cussed in the previous section. The best fit to the the- 
oretical data gives: 

(B-V)= 0.890 + 0.256 [Fe/H] + 0.048 [Fe/H]^ 
(at Mv=6.0) 

which practically overlaps the relation given by Gratton 
et al. (1997a) except for a difference in the zero point by 
A(B-V)= 0.014. As an alternative approach to the effect of 
metallicity one finds for the MS magnitude at (B-V)=0.6: 

Mv= 3.889 - 2.160[Fe/H] - 0.509[Fe/H]2 
(at (B-V)= 0.6) 

One can use either of these two relations for correcting the 
metallicity effects in data in Fig. 3, by shifting the observa- 
tional data to the [Fe/H] value of the computed isochrone. 
However, one finds that such a procedure does not improve 
the fitting given in Fig. ^, as an evidence that parallax 
errors (i.e., errors in absolute magnitudes) overcome the 
effect of metallicity in the chosen dwarf samples. 

When moving from the MS to more advanced phases. 
Fig. H gives TO V-magnitudes as a function of age for the 
various adopted metallicities. One finds that for each given 
metallicity the TO magnitudes can be arranged as a linear 
function of the logarithm of the age, with a maximum er- 
ror (for Z=0.001) not exceeding 0.4 Gyr. Table 1 gives the 
coefficients of the linear relations connecting the TO visual 
magnitudes with cluster ages (in Gyr) for isochrones with 
element diffusion taken into account. Figure shows the 
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Fig. 4. Present predictions of B-V colors at Mv=6.0 for the 
various metallicities compared with observational data given 
by Gratton et al. (1997a) on the basis of Hipparcos observa- 
tions. 

comparison between current predictions for the TO mag- 
nitudes at t=10 Gyr and for the four adopted metallicities 
with similar results in the literature. As already predicted 
in paper I, one finds that present evaluations tend to de- 
crease the cluster age for any given TO magnitude, the 
differences in Fig. |^ being mainly the consequence of the 
difference already discussed in Paper I concerning theo- 
retical luminosities. In passing we note that at the larger 
metallicities present magnitudes do not overlap Mazzitelli 
et al. (1995) predictions, as occurred for luminosities. This 
is mainly due to a small difference in adopted visual mag- 
nitude of the Sun: MDC adopted Mvq=4.79 mag. against 
our preferred canonical value Mv0=4.82 mag. 



Fig. 5. Visual magnitude of the TO isochrones as a function of 
the logarithm of the age, for the labeled metallicities for models 
with (solid line) and without (dashed line) element diffusion. 



Fig. 6. TO visual magnitudes for t=10 Gyr as a function of 
the metallicity for present models as compared with similar re- 
sults available in the current literature. For DCM 1997 models 
CM indicates the adoption by the authors of the Canuto & 
Mazzitelli (1991) treatment of overadiabatic convection while 
MLT indicates the adoption of the usual mixing length theory. 

As already known, the TO magnitude, which is de- 
fined as the magnitude of the bluest point of the isochrone, 
is a difficult observational parameter in particular in the 
most metal poor clusters. Moreover, small variations in 
the shape of the theoretical isochrone can produce rele- 
vant variation in the magnitude of the nominal TO. This 
explains the small departures from linearity disclosed by 
data plotted in Fig. |^. In a recent paper Chaboyer et 
al. (1996) suggested a different calibrator of the globu- 
lar cluster age, as given by the visual magnitude of the 
so called BTO point, that is the point brighter than the 
TO and redder by A(B-V)=0.05. The main advantage of 



Mv(BTO) is that this point is not in the vertical turn-off 
region and thus it is easier found in the observed color 
magnitude diagram. Figure ^ shows that the visual mag- 
nitude of theoretical BTO appears much more regularly 
depending on the cluster age than Mv(TO) does, further 
supporting the use of such a parameter. 



Fig. 7. Visual magnitude at the isochrone BTO for present 
models with diffusion as a function of the logarithm of the age, 
for the labeled metallicities. 



4. Horizontal branches. 

As a preliminary point, let us recall that in Paper I 
HB models have been computed assuming initial He core 
masses from a 0.8 Mq progenitor. Strictly speaking, ac- 
cording to data given in the same paper, this implies an 
age of the order of 11 Gyr for the most metal poor clus- 
ter, increasing up to 16 Gyr for Z=0.006. In Paper I we 
already discussed this issue, reaching the conclusion that 
the expected variation of HB luminosity can be neglected 
over a reasonable range of ages. Now we add that detailed 
numerical experiments show that a variation of the clus- 
ter ages of ± 5 Gyr gives a variation in the predicted HB 
magnitude of ± 0.02 mag. According to the calibration 
given in Fig. ^ a similar error in the cluster distance mod- 
ulus would produce an error in the cluster age not larger 
than 0.3 Gyr which, however, could be easily taken into 
account whenever a better precision would be required. 

Table 2 gives details for HB structures in the obser- 
vational plane at the various metallicities. Selected evolu- 
tionary quantities for the new HB models with Z=0.0006 
are given in the Appendix of this paper. A comparison 
of theoretical predictions for HB with Hipparcos data has 
been already presented in Paper I. Figure ^ compares the 
magnitude of the ZAHB (at logTe=3.85) with previous 
evaluations on the matter. One recognizes the not negligi- 
ble increase in the predicted HB luminosity induced by the 
improved physics as well as the fair agreement, toward the 
lower metallicities, with the recent computations by Caloi 
et al. (1997). As already predicted (Castellani et al. 1991) 
one finds that the dependence of the ZAHB luminosity on 
the metallicity increases with metallicity. The best fit of 
data in Fig. || gives: 

Mv = 0.993 + 0.461 [Fe/H] + 0.087[Fc/H]2 

which reproduces the theoretical predictions for ZAHB 
with diffusion within less than 0.01 mag all over the as- 
sumed range of metallicities. If diffusion is not taken into 
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account the above magnitude has to be decreased by about 
0.04 mag. 



as age indicator for galactic globulars. This is shown in 
Fig. where we report AV(TO-HB) as a function of the 
age for the four selected assumptions about stellar metal- 
licity. 



Fig. 8. Visual magnitude of the ZAHB at logTe=3.85 as a 
function of the metallicity for present models with (solid line) 
and without diffusion (dashed line) compared with values for 
the same quantity available in the literature. CCP91 indicates 
models by Castellani, ChiefB & Pulone (1991). 

For metallicities lower than Z=0.001 the above rela- 
tion can be approximated by a linear relation, as usually 
adopted in the literature. We find: 
Mv= 0.18[Fe/H] + 0.74 (no diffusion) 

Mv= 0.18[Fe/H] + 0.77 (diffusion). 

These results can be usefully compared with the fairly 
large amount of observational relations presented in the 
literature: 



Mv= 


0.15[Fe/H] - 


h 1.01 


(Carney et al. 1992) 


Mv= 


0.15[Fc/H] - 


F 0.73 


(Walker 1992) 


Mv= 


0.15[Fc/H] - 


h 0.84 


(De Santis 1996) 


Mv= 


0.18[Fe/H] - 


h 0.74 


(Gratton et al. 1998b) 


Mv= 


0.19[Fc/H] H 


h 0.97 


(Clementini et al. 1995) 


Mv= 


0.23[Fe/H] - 


h 0.83 


(Chaboyer et al. 1998) 


Mv= 


0.30[Fe/H] - 


1- 0.94 


(Sandage 1993) 



It appears that present predictions show a dependence on 
the metallicity in reasonable agreement with the evalu- 
ation by Walker (1992) and the more recent evaluation 
given by Gratton et al. (1998b) based on Hipparcos re- 
sults, and definitely smaller than required by Sandage in 
his scenario for explaining the Oosterhoff dichotomy. How- 
ever, before comparing the zero point of the magnitudes 
one has to recall that our previous relations refer to the 
ZAHB luminosity level, whereas observational data refer 
to the mean luminosity of the HB at the color of the RR 
Lyrae Gap. The connection between the two luminosity 
level has been already discussed in several paper (see Ca- 
puto et al. 1987, Carney et al. 1992, Cassisi & Salaris 
1997). One can safely assume AMv=0.08 mag as a suit- 
able estimate of the difference in magnitude between RR 
and ZAHB. Thus our zero points will become Mv=0.69 
(diffusion) or Mv=0.66 (no diffusion), respectively. One 
finds that we are predicting HB with the same depen- 
dence on metal content and only 0.05 mag. brighter with 
respect to the recent estimates by Gratton et al. (1997). 

One can finally connect theoretical results concerning 
HB stars with previous predictions about the TO magni- 
tudes to give a theoretical calibration of the difference in 
magnitude between HB and TO, AV(TO-HB), often used 



Fig. 9. The difference in magnitude between HB and TO as 
a function of the age for present models with (solid line) and 
without (dashed line) element diffusion for the labeled values 
of metallicities. 



5. The role of the original He content. 

All the theoretical computations referred in the previous 
sections assume an amount of original He content equal 
to Y=0.23, which is, at the present, the most popular and 
widely adopted estimate for the original He abundance in 
metal poor. Population H stars. However, in Paper I we 
already discussed the large uncertainty in the evolution- 
ary determination of this parameter. According to such 
an evidence, we will investigate the role played by the as- 
sumption about Y on the present theoretical scenario. To 
this purpose we reinvestigated the evolutionary behavior 
of the Z=0.001 models (with element diffusion) but under 
the two alternative assumptions Y=0.21 or Y=0.25. 



Fig. 10. Present isochrones in the CM diagram for Z=0.001, 
and the labeled values of age and original He abundances. 

Figure |l^ compares selected isochrones for Z=0.001 
and for the three adopted values of helium abundance. As 
a result one finds the following relations: 

unevolved MS (at B-V=0.6): 
AMv/AY--2.5 (slowly depending on (B-V) 

TO (at t=10 Gyr): 

AMv= 0.019 (Y=0.21^0.23) 

AMv=0.034 (Y=0. 23-^0. 25) 

HB (at logTe=3.85): 

AMv= -0.084 (Y=0.21-=-0.23) 

AMv= -0.068 (Y=0.23~0.25) 

As a whole, present results (where element sedimenta- 
tion is taken into account) appear in reasonable concor- 
dance with canonical evaluations given by Rcnzini (1991) 
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for the TO magnitudes and by Buzzoni et al. (1983) or, 
more recently, by Bono et al. (1995) for HB luminosities. 
One finds that a variation of Y within the assumed limits 
(Y=0.23±0.02) gives a maximum variation of 0.034 mag. 
in the TO magnitude and of 0.10 mag. in the difference of 
magnitude between TO and HB. According to the calibra- 
tion given in previous sections this uncertainty implies an 
error of about 0.2 Gyr for ages from the TO magnitude, 
and an error of about 1 Gyr for ages from the "vertical" 
method, i.e., from the difference in magnitude between 
TO and HB. 

From the previous sections one also derives 
AMv(TO)- 0.37 A[Fe/H] and AMv=0.18 A[Fe/H]. As- 
suming an uncertainty of ± 0.2 dex on current evaluations 
of cluster metallicity, we obtain that such an error drives, 
by itself, a variation of, about, 0.08 mag in the predicted 
TO, and of about 0.04 mag on the difference in magnitude 
between TO and HB. As a result, one finds that with the 
assumed uncertainty in both Helium and metal content, 
even with perfect photometry, the age of a cluster cannot 
be determined better than ± about 0.8 Gyr from the TO 
magnitude or better than, about, 1.4 Gyr from the vertical 
method. 

We remind that in all cases we assumed scaled solar 
composition as given by Grevesse & Noels (1993). How- 
ever, for metallicity not too much high, the effect of the 
a-enhancement on the evolutionary tracks and isochrones 
can be simulated by using a scaled solar mixture of to- 
tal metallicity equal to the actual one (see Salaris et al. 
1993). Detailed calculations recently performed by Salaris 
& Weiss (1998b) show that the a-enhancement does not 
influence at all the AV(TO-HB) parameter for metallic- 
ities up to ZwO.004, while for Z=0.01 the variation of 
AV(TO-HB) is about 0.1 mag. Thus the influence of the 
a-enhancement on AV(TO-HB) is expected to be negli- 
gible up to our highest adopted metaUicity (Z=0.006). In 
the case of a-enhancement we only expect, for Z=0.006, a 
shift of the color of the RG branch, well within the known 
uncertainties due to the uncertainty on the efficiency of 
superadiabatic convection and on the adopted models at- 
mospheres. 



Fig. 12. Isochrones for age between 9 and 11 Gyr and ZAHB 
compared to the CMD of M5 (data from Sandquist et al. 1996). 
Composition, distance modulus and reddening as labeled. The 
adopted mixing length is 2.3 Hp. 



MS at a given color is in principle an indicator of the 
original He content independent on the cluster age and 
not affected, as the well known parameter R is, by the 
still large inaccuracy on the rate of the ^^C-a reactions. 
This indicator is becoming more and more relevant vis-a- 
vis the increased capability of precise photometry of faint 
MS stars in galactic globulars. Combining the previous re- 
lations one finds AMv(HB-MS)-6.5AY. Thus AMv(HB- 
MS) appears rather sensitive to the amount of original He, 
increasing by 0.07 mag. when Y is increased by only 0.01. 
However one should remember that significant theoretical 
uncertainties, due to the uncertainty in the efficiency of 
superadiabatic convection and to the choice of color trans- 
formations, together with observational uncertainties, due 
to the evaluation of the cluster reddening and metallicity, 
make de facto this quantity a difficult parameter to be 
used as helium calibrator. 



Fig. 13. Theoretical predictions concerning the HB magni- 
tudes from the present paper (full line: diffusion) or from CCP 
(dashed line) compared with recent observational estimates as 
derived either from the subdwarf fitting (open symbols) or from 
RR Lyrae parallaxes (filled symbols). 



6. Discussion and final remarks 



Fig. 11. Isochrones for ages between 10 and 12 Gyr and ZAHB 
compared to the CMD of M68 (data from Walker 1994). Com- 
position, distance modulus and reddening used for the fit as 
labeled. The adopted mixing length is 2.0 Hp. 



As already discussed by Caputo et al. (1983) one fi- 
nally finds that MS and HB have opposite behavior as far 
as the amount of original He is concerned. This implies 
that the difference in magnitude between the HB and the 



Before closing the paper let us use the presented theo- 
retical scenario to fit the C-M diagram of selected glob- 
ular clusters. For a first test we chose C-M diagram pre- 
sented by Walker (1994) for the metal poor cluster M68. 
A preliminary discussion on that matter has been recently 
presented by Brocato et al. (1997) on the basis of Paper 
I evolutionary results and using Kurucz (1992) transfor- 
mations. Without repeating the discussion given in that 
paper we present in Fig. 11 the best fitting of the observa- 
tional diagram as obtained on the basis of our Z=0.0002 
isochrones with element diffusion, adopting the very last 
version of Kurucz's model atmosphere (Castelh 1997a, b), 
where the new solar abundances and the enhancement of 
a-elements have been taken into account. 
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The best fitting is achieved for an age of about 11 
Gyr, with distance modulus and reddening as labeled. 
Note that to fit the color of the RG branch we adopted 1= 
2.0Hp. However, Fig. |ll| shows that with such an assump- 
tion the MS appears slight bluer than observed, requiring 
lower values of the mixing length parameter and indicating 
that the CM diagram of the cluster can hardly be fitted 
with precision assuming a constant mixing length for all 
the evolutionary phases. By neglecting this (false) prob- 
lem of temperature of cool stars, one finds that the best 
fitting requires a cluster distance modulus DM= 15.30, in 
excellent agreement with Gratton et al. (1997) who used 
Hipparcos parallaxes to find for the cluster DM=15.31 and 
a (mean) age of 11.3 Gyr. 

As a test of theory at larger metallicity. Fig. ^ 
present the best fitting of the C-M diagram presented 
by Sandquist et al. (1996) for the intermediate metal- 
licity cluster M5. Now the best fitting gives an ages of 
about 10 Gyr and a distance modulus DM= 14.54. On 
the basis of Hipparcos parallaxes Gratton et al. (1997) 
give DM=14.60±0.07 with a (mean) age of 10.5 Gyr. 
On the same ground, Chaboyer et al. (1998) gives DM= 
14.51±0.09 with an age of 8.9±1.1 Gyr. 

The presented evolutionary models appears in both 
cases to give excellent agreement with independent evalu- 
ations of the cluster distance moduli based on Hipparcos 
data. This strongly increases the confidence in HE stars 
as (theoretical) standard candles and, in the same time, in 
the reliability of the derived cluster ages. One may notice 
that both clusters we are dealing with have been recently 
fitted with less updated physics, corresponding to step 4 
in Paper I. In that case it was derived t= 12.2 Gyr for 
M68 (Safaris et al. 1997) and, with slight different assump- 
tions about the cluster chemical composition (Y=0.235, 
Z=0.0015), t= 10.9 Gyr for M5 (Salaris & Weiss 1998a). 
Comparison with present results casts new light on the 
further rejuvenation of cluster ages induced by both the 
subsequent updating of the physics and the introduction 
of the element diffusion. 

One can safely assume a conservative error of ± 0.1 
mag in our estimates of the difference in magnitude be- 
tween TO and HB, due to the arbitrarity introduced when 
the theoretical HB and isochrones are fitted to an observa- 
tional CM diagram. As an example, we note that Brocato 
et al. 1997 using the theoretical models of paper I, find 
for M68 a distance modulus (DM=15.25) slightly different 
from the present result, due to the slightly different way in 
which the theoretical HB is fitted to the observational one. 
This variation in the distance modulus together with the 
adopted color transformations (which influence the look 
of the fit) yields a difference of ~ 1 Gyr in the estimated 
age. With the above quoted assumptions about the uncer- 
tainty in the chemical composition we conclude that the 
adopted theoretical scenario gives for our clusters: 



M68: t= 11 ±1.4 ±1.0 Gyr 
M5: t= 10 ±1.4 ±1.0 Gyr 

where the first error is due to the uncertainty in the chem- 
ical composition, while the second represents the uncer- 
tainty in the fit. In passing, we note that the close sim- 
ilarity between the C-M diagrams of M68 and M92 (as, 
e.g., discussed in Brocato et al. 1997 and Salaris et al. 
1997) drives to the conclusion that both clusters should 
have quite similar ages. According to the above discussion 
we suggest for these very metal poor globulars an age of 
the order of 11 Gyr against the 14 Gyr recently derived 
by Pont et al. (1998). One may notice that the above age 
evaluations could suggest a possible correlation between 
cluster ages and metallicity, the more metallic cluster be- 
ing also the younger one. The evidence from the figure 
appears clear enough, however the source of possible er- 
rors in the photometry and in the fitting do not allow firm 
conclusions about a problem which deserves much more 
accurate investigations. 

As a final remark, let us here remind once more that 
the above age estimates rely on the theoretically predicted 
HB luminosity. We have already quoted the good agree- 
ment of such a prediction with cluster distance moduli as 
derived by the fitting of Hipparcos subdwarf magnitudes. 
However one has also to remind that several estimates 
of RR Lyrae luminosities based on Hipparcos data give 
sensitively fainter magnitudes. As shown in Fig.|l3|, the is- 
sue is far from being clearly settled. Here we can only say 
that if such faint magnitudes will be eventually confirmed, 
present theory is obviously overestimating the He cores at 
the end of RG evolution, likely as a consequence of a cor- 
responding overestimate of the efficiency of cooling along 
the RG phase. In this case, data in Table 1 indicate that 
cluster ages should be increased by Alogt^O.4 SMv, where 
SMv represents the difference between the actual and the 
predicted magnitudes of RR Lyrae stars. 
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8. Appendix 

Table 3 lists selected evolutionary quantities for the new 
models with Z=0.0006 when the efhciency of element sed- 
imentation is neglected. Left to right one finds: mass of 
the model, age, luminosity and effective temperature at 
the track Turn-OfF (TO), mean luminosity of the RGB 
"bump", age at the He flash and the corresponding lu- 
minosity, mass of the He core and amount of extra-He 
brought to the surface at this time. Table 4 gives the same 
quantities but for models where element sedimentation is 
taken into account according to the procedure reported in 
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Paper I. Table 5 gives details on the isochrone TO over the 
explored range of ages for the new models with Z=0.0006. 
Table 6 gives selected evolutionary quantities for the new 
HB models with Z=0.0006. 

Finally we compare the evolutionary characteristics of our 
models with those of similar models, for the same chem- 
ical composition and slightly different physics, recently 
presented by Straniero et al. (1997) finding a substan- 
tial agreement. Regarding the H burning phases, there is 
a good agreement in the luminosity and effective temper- 
ature of the TO region, while our TO ages are lower by 
about 10%, presumably due to the adoption of a different 
EOS. As far as the evolutionary characteristics at the He 
flash is concerned, the value of the extrahelium brought to 
the surface during the first dredge-up is the same, while 
our He cores are larger of about 0.005 M0 . This very slight 
difference could be due to our adoption of updated calcula- 
tions for plasma neutrino energy losses (Haft et al. 1994). 
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Table 1. Coefficients of tfie linear relation: Mv"^*^ = aLog(age) + b which connects the visual magnitude at the isochrone TO 
(Mv"^") with the TO age for present models with element diffusion and Y=0.23. Age is in Gyr, metallicity as indicated. 



z 


a 


b 


0.0002 


2.70 


0.98 


0.0006 


2.55 


1.30 


0.001 


2.48 


1.46 


0.006 


2.33 


1.83 



oo 
On 

> 
O 

o 

(N 

^ \ M 0.530 0.540 0.550 0.560 0.580 0.600 0.620 0.650 0.670 0.690 0.700 0.720 0.750 0.770 0.800 

^ . Mv- 3.954 3.658 3.415 - 2.806 2.436 2.045 1.457 - - 0.663 0.499 0.423 0.427 0.439 

■ B-V -0.214 -0.196 -0.183 - -0.150 -0.132 -0.112 -0.081 - - 0.008 0.076 0.233 0.327 0.446 
CN Z=0.0006 YijB=0.229 
m Mc= 0.509 Me 



Table 2. Selected evolutionary quantities in the observational plane for ZAHB models with Y=0.23 and element diffusion. 
Masses are in solar units. 



Z=0.0002 Yhb=0.226 
Mc= 0.515 Me 



M 0.530 0.540 0.550 0.560 0.580 0.600 0.620 0.650 0.670 0.690 0.700 0.720 0.750 0.770 0.800 

yVv 3.762 - 3.084 - - 1.948 - 0.650 0.525 0.516 - - 0.523 - 0.487 

B-V -0.202 - -0.162 - - -0.108 - 0.027 0.139 0.331 - - 0.612 - 0.660 

^ ■ Z=0.001 Y//B=0.229 



OO 
ON 



M<,= 0.507 M 







o ■ 

' M 0.530 0.540 0.550 0.560 0.580 0.600 0.620 0.650 0.670 0.690 0.700 0.720 0.750 0.770 0.800 

^ ; 3.673 - 3.099 - - 1.702 - 0.548 - - 0.590 - 0.555 

C3 B-V -0.202 - -0.171 - - -0.102 - 0.187 - - 0.643 - 0.718 



Z=0.006 Yhs=0.231 
Mc= 0.500 Me 



H ; M 0.530 0.540 0.550 0.560 0.580 0.600 0.620 0.650 0.670 0.690 0.700 0.720 0.750 0.770 0.800 
5r \ M_v 3.101 2.748 2.368 1.904 0.793 0.833 _ _ _ _ 0.781 - - - 0.747 

B-V -0.157 -0.138 -0.118 -0.094 0.360 0.610 _ _ _ _ 0.898 - - - 0.942 



2 



Table 3. Selected evolutionary quantities (see text) for stellar models with Y=0.23 and Z=0.0006 without element diffusion. 
Ages are in Gyr and masses are in solar units. 



M 




LogL^° 


LogTe^° 


LogL 


^flash 


LogL^'"'''' 


Mc 


J- HB 


0.6 


34.9 


-0.083 


3.762 




38.56 


3.375 


0.5073 


0.232 


0.7 


19.4 


0.107 


3.785 


1.956 


22.40 


3.389 


0.5083 


0.237 


0.8 


11.5 


0.278 


3.808 


2.019 


13.85 


3.384 


0.5062 


0.242 


0.9 


7.4 


0.460 


3.834 


2.091 


9.122 


3.365 


0.5019 


0.2455 


1.0 


5.1 


0.638 


3.870 


2.133 


6.33 


3.372 


0.5022 


0.249 



Table 4. The same as in Table M for models with element diffusion. 



M t'^° LogL^Q LogTe^° LogL'"""' t-^'""'^ LogL^'""'' Mc Yhb 



0.6 


31.9 


-0.137 


3.752 




36.10 


3.376 


0.5192 


0.200 


0.7 


18.1 


0.071 


3.777 


1.917 


21.27 


3.391 


0.5142 


0.217 


0.8 


10.8 


0.251 


3.802 


1.990 


13.36 


3.388 


0.5094 


0.229 


0.9 


6.95 


0.428 


3.829 


2.065 


8.88 


3.382 


0.5063 


0.237 


1.0 


4.9 


0.619 


3.862 


2.126 


6.20 


3.374 


0.5037 


0.243 
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Table 5. Selected TO quantities for H burning isochrones with Y=0.23 and Z=0.0006, with and without element diflFusion for 
the labeled values of age (Gyr). 



t™ 


LogL^o 


LogTe^° 


M"^"^ [Mq] 


standard 


8 


0.510 


3.831 


0.893 


9 


0.461 


3.824 


0.865 


10 


0.418 


3.818 


0.840 


11 


0.373 


3.812 


0.818 


12 


0.320 


3.807 


0.798 


13 


0.288 


3.803 


0.781 


14 


0.264 


3.800 


0.766 


15 


0.242 


3.797 


0.753 


16 


0.216 


3.794 


0.740 


17 


0.192 


3.791 


0.728 


18 


0.169 


3.789 


0.717 


diffusion 


8 


0.473 


3.823 


0.883 


9 


0.420 


3.816 


0.854 


10 


0.381 


3.809 


0.831 


11 


0.339 


3.803 


0.809 


12 


0.294 


3.799 


0.789 


13 


0.265 


3.795 


0.772 


14 


0.238 


3.791 


0.757 


15 


0.207 


3.787 


0.743 


16 


0.184 


3.784 


0.730 


17 


0.162 


3.781 


0.718 


18 


0.122 


3.778 


0.706 
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Table 6. Selected evolutionary quantities for ZAHB models with Z=0.0006. 



standard 



Z=0.0006 Yhb=0.242 Mc= 0.5062 M, 



M [Mo] 


LogL 


LogTe 


tife [Myr] 


0.53 


1.347 


4.355 


95.8 


0.55 


1.392 


4.275 


91.4 


0.60 


1.548 


4.127 


85.3 


0.65 


1.681 


3.938 


81.5 


0.70 


1.743 


3.770 


79.3 


0.75 


1.786 


3.731 


78.8 






diffusion 




Z=0.0006 Yhb 


=0.229 Mc= 


: 0.5094 Mo 


M [Mo] 


LogL 


LogTe 


tife [Myr] 


0.53 


1.353 


4.364 


94.5 


0.55 


1.394 


4.276 


90.2 


0.60 


1.526 


4.131 


84.4 


0.65 


1.663 


3.938 


80.9 


0.70 


1.726 


3.775 


79.4 


0.75 


1.762 


3.734 


76.8 


0.80 


1.786 


3.724 


76.7 



